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ABSTRACT 

The claim of the inflation theory having explained large scale flatness and 
absence of monopoles and strings is examined from the viewpoint of the observed 
scales having originated from very small ones, on which the density fluctuations 
of the curvaton and relics are inevitably of order unity or larger. By analyzing the 
density evolution in two different gauges to ensure consistency of the smoothest 
possible pre-inflationary component, viz. radiation, it is found that the 0(1) 
thermal fluctuations on the thermal wavelength scale presents a problem to the 
linear growth theory (to the best of the author's knowledge this is also the first 
time the longitudinal and synchronous gauge are formally reconciled in detail by 
explicit gauge transformation). Even if one invokes a classical thermodynamic 
argument to ensure the start of inflation and the significance of the role of relic 
radiation by the time such 0(1) fluctuations exit the horizon, difficulties may 
still resurface in the coherent oscillation stage of reheating. 



1. Introduction 

The theory of inflation (Guth (1981); Albrecht & Steinhardt (1982); Linde (1982)) 
is thought to have provided attractive solutions to a number of cosmological problems, in- 
cluding the large-scale homogeneity and flatness of the universe. It garnered strong support 
from the COBE and WMAP observations of the cosmic microwave background radiation (CMBR) 
(Smoot et al (1992); Bennett et al (2003); Spergel et al (2007); Hinshaw et al (2009); 
Komatsu et al (2011)). Quantum fluctuations in a scalar inflaton fleld can explain the origin 
and near scale-invariant spectrum of the primordial density perturbations, although getting 
the amplitude right requires flne tuning. 

Apart from the central question of structure formation, however, we wish to revisit some 
of the other key issues of cosmology that inflation is claimed to have satisfactorily addressed, 
as they too lend as crucial evidence in favor of the theory if the logic is sound. In particular. 
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the reason for the 'global' flatness of space and the absence of strings and monopoles is, 
in the context of inflation, due to the ability of the rapid expansion in drastically diluting 
these features. This means, at the onset of inflation the density of curvatons and relics 
was probably on par with that of the inflaton. On the relevant cosmological scales that 
were deeply embedded within the horizon at the time, the fluctuations in the densities of 
the unwanted components are then initially very large, and are likely to be governed by 
quantum effects. In general, it is reasonable to assume that the smoothest component is the 
one with the largest number density, viz. the massless particles. Therefore in this paper we 
examine closely the state of a pre-existing photon population on small scales, by taking the 
universe before inflation as radiation dominated. This will set lower limits on the severity of 
the initial fluctuations. 

Indeed, early versions of inflation theory did envisage a pre-inflationary Priedmann- 
Robertson- Walker (FRW) stage. One of the claimed advantages of the theory was that 
inflation effectively erased all traces of that 'pre-historic' phase, although it was realized early 
(Frieman & Turner (1984)) that this is not strictly true; what it does is not to eliminate 
perturbations but to stretch them to unobservable scales. This leaves open the question of 
whether there could be perturbations on very tiny scales that are stretched to observable 
size. It has been shown (Magueijo & Singh (2007)) that without drastic modifications such 
perturbations could not explain the power spectrum of perturbations. Here we wish to argue 
that under rather general conditions these perturbations would both be inevitable and far 
too large to be consistent with either perturbation theory or observation. This is because 
thermal radiation has 0(1) fluctuations on scales comparable with the thermal wavelength, 
and these scales have been stretched to cosmologically relevant ones or beyond. 

In yet another manner of speaking, while it has previously been argued (Vachaspati & 
Trodden (1998)) that inflation requires sufficiently fine-tuned initial conditions to ensure 
ultra-smoothness in the relevant part of the universe at the onset of the process, we suggest 
that it is essentially impossible to satisfy that condition if there is a pre-inflationary FRW 
stage consisting primarily of radiation or other even less homogeneous components. 

To focus the readers' attention upon the effect we wish to address, we shall here-and- 
after consider inflation driven by a scalar field of small density fluctuations Sp.,, relative to 
the radiation Spr, and with an equation of state close to p„ = — Our conclusion is actually 
insensitive to such assumptions. We note in passing that fluctuations of an inflaton with 
vanishing 1 -\- Wy (or 0^) do not necessarily lead to infinite (or indeterminate) curvature 
perturbations ( as each scale leaves the horizon, i.e. this claim is true only if there are no 
other accompanying components present. 
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2. The model 

We assume that in the pre-inflationary phase the segment of the universe from which 
today's Hubble volume evolved contained a fluid of massless particles, with interactions weak 
enough to treat them as free, but strong enough to maintain thermal equilibrium. We also 
assume the existence of an inflaton field providing vacuum energy that would eventually 
come to dominate and induce a period of inflation. 

The background space-time is taken to be homogeneous, isotropic and spatially flat. 
We consider only scalar metric perturbations, flrst in the longitudinal or Newtonian gauge, 
then in the synchronous gauge, then cross compare the results to flnd consistency. In the 
former, we assume that the anisotropic stress is negligible, i.e. the two invariant potentials 
are equal: $ = ^. Thus we can write the perturbed metric as 

ds^ = a\T)[{l + 2^)dT^ - (1 - 2^)dx'dx'], (1) 

where r is the conformal time {adr = dt). In the latter the metric is 

ds^ = a^(T) [dr^ - (5ij + hij)dx'dx^] , (2) 

where 

hij = -ihSij - 6{didj - \5ijV'^)r] (3) 

are respectively the trace and traceless contributions to hij (Ma & Bertschinger (1995)). 

We suppose that the vacuum energy starts to dominate at an initial time t\ at which 
Pr = Pv, where stands for the energy density of the radiation and pv for that of the 
inflaton fleld. Inflation ends at a reheating time Trh when the vacuum energy is converted 
to radiation. Since Pv is nearly constant during the inflationary era, this means that the 
physical temperature Tphys of the radiation is approximately the same at ri as the reheating 
temperature just after Trh- This is a large temperature, but still some orders of magnitude 
below the Planck energy. It is convenient to use a 'comoving temperature' T — aTphys which 
is nearly constant outside the inflationary period. More precisely, g^^^T is constant, where 
g is the number of helicity states of massless particles (with fermions counted as |); for the 
universe since reheating, allowing for the effect of neutrino decoupling and electron-positron 
annihilation, 

= (l^y Vo, (4) 

where Tq — 2.7K is the CMB temperature today. Thus the approximate equality of the 
physical temperatures before and after inflation means that 

^rh 7'(7"rh) «rh a(7"rh) 
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where N is the number of e-folds of inflation. Equivalently, the comoving temperature in 
the pre-inflationary phase is 

Ti « O.Se-^To. (6) 

It is now necessary to be precise about the power spectrum of thermal radiation. On 
scales above one thermal wavelength, i.e. k <T it is given by standard thermodynamics as 

n^) - - (7) 

In particular, V{k) is 0(1) on the wavelength scale k ^ T (note also that since both k and 
T are 'comoving', k/T = Ajphyg/Tphys)- This is important because k ^ T length scales exited 
the horizon at an early stage during inflation. The situation for the even smaller scales of 
k > T that exited the horizon later is more complicated and requires a separate treatment, 
to be given in the following section. 



3. Thermal perturbations on sub-wavelength scales 

Owing to their immense importance we calculate from flrst principles the density per- 
turbations in a system of massless particles in thermal equilibrium at temperature T. Of 
course, in the classical region where /c 3> T, the answer is well known and easily derivable 
from thermodynamic arguments (see e.g. Ferreira (2008)). But the result for k <T is not 
so obvious and not so simple to interpret. 

For simplicity we consider only a single massless scalar field 0. For the general case, 
we simply have to multiply by the number of helicity states g. It is sufficient to consider 
Minkowski space. In equilibrium at temperature T, the mean occupation number for the 
mode with wave number k is n(k) = {e^^'^ — 1)"^. 

The energy density of the field is u{t, r) = | : [0^(t, r) + (V<;6)^(t, r)] :, where the colons 
denote normal ordering. In terms of creation and annihilation operators, 

- a*(k)a*(k')e-'(''+'''>"- + 2a*(k)a(k')e^('^-'^'>"-] , (8) 
For the mean value of u only the last term in (8) contributes, so 

^""^ ^ (2^ y ^'^^^(^) = ^ / kMk)dk = -T\ (9) 
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Of course the result for photons is the same, except for a factor 2 arising from the two 
polarization states. 

What we want to calculate is the dimensionless density perturbation power spectrum, 
V{k) = ^^P{k) = d'r^{v)e-^^-, (10) 

where ^, the density correlation function, is defined by 

«r)^ <"'°-f"f°'> -l. (11) 

To evaluate this, we now substitute (9) into (w(0, r)u{0, 0)) twice. Clearly, the only terms 
that will give nonzero contributions are those with two as and two a*s. Moreover, since there 
is no correlation between different wave vectors, there must be two delta functions relating 
the arguments of the as to those of the a*s. One term will simply reproduce the square of (u), 
so cancels the 1 in (11). In addition, there is a divergent term that would be present even at 
zero temperature, which must also be subtracted. In the remaining terms, the integration 
variables can be chosen so that each contains the same delta function (53(k — ki — k2): 

- 2(2^ j^j ^^3(k-k.-k.).(/.0KM + l] 

X [{kik2 + ki • ks)' + {kik2 - ki • k2)']. (12) 

Next we can do the angular integrations, using ki ■ k2 = \{k'^ — k\ — A;|). It is convenient to 
introduce the two dimensionless integration variables X — {ki-\- k2)/k and Y — {ki — k2)/k. 
Thus we find 

r{k) = — TTTT? / dX dY 



2567r4(M)2 (,{X+Y)k/2T _ 1 1 _ e-(X-Y)kl2T 

X [(X2-1)2+(F2_i)2|_ ^^3^ 

To perform the integrals, we may expand the denominators, obtaining 

^{k) = 900 (-^) dX e-^'^'^'''/'^ / dy e-('-^)^'=/^^ 

X [{X^ - If + {Y^ - If]. (14) 

(Note the differing ranges of the summations. The s = term arises from the +1 in (12).) 
The integrals can now be performed. To write the result in a reasonably concise form, we 
introduce the abbreviation 

/. = ^, (15) 



-6- 



Then we find 

8 0<j 'DO 



^ e^-{f^_s - ^fr-s + 3) - e-f^-'{fl, + 3/,_3 + 3) ^^^^ 

fr+s fr 



I r—s 



We are most interested in the limiting forms of V for small and large k/T. Expanding 
in powers oi k/T we find that all terms in {k/TY cancel, and the leading term is 



W~^^1.(;^ = ^C(4) = ^ (k^% (17) 



in agreement with the classical expression. To further elaborate upon this agreement, (17) 
is smaller than the power spectrum of thermal photons in the small k limit by a factor 
of two, which is exactly what one expects because photons have two internal degrees of 
freedom [i.e. {5uY and u are both larger, hence {5u/{u)Y smaller, by two times relative to 
the massless particles). 

So far as the large-/c behaviour is concerned, we can drop all terms involving decreasing 
exponential factors. The only surviving terms are those with s — Q and where the two 
exponentials cancel. Keeping only these, we have 

The leading term is the one in 1//^, which gives 

900C(4)A;^ lOfc^ ^ 
nk)- ^sri =^ {k^oo). (19) 

It is interesting to note that the point at which the two asymptotic forms (17) and (18) agree 
is where k/T — 6. 

Thus far we have considered a static universe. Suppose however that we have an ex- 
panding universe in which expansion is slow enough to allow the radiation to remain close 
to thermal equilibrium. The results above will then hold to a good approximation. 

The physical interpretation of this result is far from obvious. It suggests that 'P{k) 
becomes much larger than unity when k ^ T. This means that the standard deviation of 
a measurement of the radiation energy in a volume smaller than the typical wavelength of 
the radiation is much larger than its mean value. Given that the energy cannot be negative. 
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this implies a highly non-gaussian distribution. Clearly it would be difficult to devise an 
experiment to measure directly the energy density perturbations on scales smaller than the 
typical wavelength of the radiation. We are here interested in a 'measurement' mediated by 
perturbations in the gravitational field. 

Here we are discussing a quantum fluctuation in the energy of the radiation gas. One 
would normally assume that there is a process of decoherence, induced by the exponential 
expansion, in which the quantum distribution is replaced is replaced by a classical statisti- 
cal distribution. In estimating the magnitude of the density perturbations as the universe 
expands, we might use a 'typical' value of S. But here the square root of V{k) cannot be 
regarded as a 'typical' value of 6 on the length scale 1/k. 



4. Evolution equations: the longitudinal gauge 

We turn to the question of how the thermal radiation perturbations on the scale of one 
thermal wavelength or larger, (17), will evolve through the inflationary era. 

To begin with, the conformal Hubble parameter T-L — a/a — aH where a = da/dr is 
given by 

-H^ = (20) 

where p = Pr+Pv is the total density. Of course, pr while in the early stages of inflation 

we may assume that Pv is a constant. It is then convenient to introduce the dimensionless 
variable 

y - - (21) 

so that 

(This equation would not hold in the late stages of inflation, when the scalar inflaton fleld (p 
starts to roll down towards its minimum, 0^ is no longer negligible compared to the potential 
energy V{(f)), so the 'vacuum' energy becomes inhomogeneous and = P/p starts to deviate 
from |, indeed becoming large and negative. For the moment we exclude that era; we will 
not need explicit solutions there, because wc can use a conservation law.) 

Now the gauge-invariant potential $ satisfles a second-order evolution equation (see for 
example Peter & Uzan (2009)), 



+ 3(1 + cl)H^ + [2H + H\l + 3c2)]$ - c^V^^ = 0. 



(23) 
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If we are interested in a universe that contains only radiation and vacuum energy, in which 
the density fluctuations in the latter may be neglected as part of our ansatz (see the end of 
section 1) we may set — \, obtaining 



s 

/2\a^ , 11,2 



$ + 47^$ + 2(?^ + 7^^)$ + = 0. (24) 

This equation too will not hold in the later stages of inflation. Denoting derivatives wrt y 
by primes, we have for any function f, f = 'Hyf- Thus (24) becomes 

n^y^^" + {bv?y + nn'y^)^' + 2{n^ + nn'y)^ + ik^^ = o. (25) 

Prom (22) we also have 

nn'y =^(y'-l-). (26) 



2 V y'y 

Substituting (22) and (26) into (25) and rearranging the terms gives 

y\y^^" + 6|/$' + 4$) + y-^{^" + Ay^') + = 0, (27) 
where we have written 

^ 37^ ^ ^^^^ 
(Note that T-Li — UiHi, but Hi is not the near-constant Hubble parameter H fa Hi^s. during 
inflation; at Ti, p — 2py and therefore H-^ — y/2HinR ~ V2H). 

In (27) we have chosen to group the terms in such a way that if $ ~ i/" they behave 
respectively like y"'~^^, y""'"^ and q'^y"'. The point at which p^. — p^ is given hj y = 1. For 
y <^ 1, we are deep in the radiation-dominated regime and the flrst bracket is small; while for 
y ^ 1 vacuum energy dominates, and the second bracket is small. The time of horizon exit 
is deflned by /c^/T^ex — 1- Since this occurs during inflation, and well after radiation- vacuum 
equality, we must have g ^ 1, and then y'^^ — 3q^. 



5. Evolution before and during inflation 

Let us flrst consider the solution well inside the radiation-dominated era, when y <^ 1. 
Then the flrst bracket in (27) is neghgible and it is straightforward to solve the equation. 
The well-known solution is 

A f sinqy \ B f cosou . \ 

2g2y2 y qy j 2q^y^ \ qy J 

where A and B are arbitrary constants. If we impose the natural condition that the pertur- 
bation remains flnite at early times, we require B — Q. 
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The density contrast is given in terms of $ by 

+ U^) + |A;2$ = -\H'^5, (30) 

or 

5r = -2[y$' + $ + g2(y2 + y-')-^$]. (31) 

Over most of the range, where l/g'<^y<^l, we can ignore terms involving inverse powers 
of qy, so (29) gives 

5r ~ Acosqy, (32) 

which represents the expected oscillation with approximately constant amplitude. The typ- 
ical value of A is given by the expected amplitude of thermal fluctuations. 

Now let us turn to the other extreme, when |/ 3> 1 and vacuum energy is dominant. 
In that case, the second bracket in (29) is neghgible. It is again possible to find an explicit 
solution to the equation, by changing variable to l/y. The solution in this case is very 
similar: 

where again C and D are arbitrary constants. 

We still need to connect these two solutions (29) and (33), to relate the constants C and 
D to A. In the intervening region, near the crossover point ?/ = 1, it is not possible to find 
an explicit solution. But we can get a good idea of the relation by simply assuming that (29) 
applies throughout the whole of the radiation-dominated era, y < 1, and (33) throughout 
1/ > 1, and matching the two sX y = 1. Assuming that q is large enough for the first terms 
in all the brackets to be negligible, the boundary conditions on $ and give us 

Acosq — C cos g — D sing, 
A sing = — C sing — D cos g, 

yielding the solution 

C = ^cos2g, L> = -^sin2g. (34) 

These relations might be modified by a more accurate treatment, but C and D would still be 
related to A and 5 by a symbolic rotation; the main effect would be to change the effective 
rotation angle from 2g. Note also that even if we did not impose the initial condition of 
B — but instead took it be comparable to A, C and D would still be of a similar order of 
magnitude to A. Prom the remarks made after (32) therefore, one expects 



$ = ^ ( ^^'^ 'i/y _ cos / I 

2gV \ (i/y ) 



(35) 
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where \5th\ denotes the thermal oscillation amplitude at the scale of interest. 

The behavior of the two terms in (33) when y > q and the perturbation is outside the 
horizon is very different. The leading terms give 

The corresponding terms in 6^ are given by (31), and are 

It is interesting that the terms in y"^ and y~^ cancel. Thus although the dominant term in 
$ is the D term, in 5 it is the C term. Equivalently the two solutions are, in terms of the 
radiation density contrast and peculiar velocity divergence 5r and = • v^, 

where 9^ is a special case of the generic equation for the total divergence 0: 

A;2($ + H^) = (7^2 - 'Ha)e = ^TiGa^ip + p)e (39) 
in the limit of = —p^ andpr = Pr/3 when {p+p)9 = '^i{pi+Pi)Oi = {Pr+Pr)&r = 4:pr9r/S. 



6. Reheating and subsequent evolution 

It is useful to define the quantity C, known as the total curvature perturbation, that has 
showed itself to be particularly useful after a mode exited the horizon: 

C = * z -^^-3H— ^ = $H ^ , 40 

3(1 + ^;) pr + Pv n'^-n ^ ^ 

Note that ( is infinite only in the case of single field inflation of flnite 6^ but vanishing 1 + Wy, 
i.e. C is finite once there are other, 1 + > components present. In terms of the variable 
y, (40) can be written as 



C = ^ + ^(y*' + $ + ^$). 



(41) 



Using the evolution equation, one can then show that 



C = -^!^. (42) 
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Both the numerator and denominator of C are small, but the quotient need not be, as we 
shall find out. 

Let us evaluate these quantities for the solutions given above in the infiationary era. 
We are interested in the values reached well outside the horizon, so we may use the late-time 
expression (36). It is easy to see that the C term yields a finite and constant hmit for C, but 
the D term is proportional to y and therefore fades away with time. More precisely 

C = -\C ^ -\Acos2q, or C = — 7= D = ~Asm2q, (43) 

with the reminder that C and D are both of order the thermal amplitude, (35). 

How does one connect with the radiation era after reheating? Here $ is given by (29) 
with A and B replaced by A and B to distinguish the amplitudes from the pre-inflationary 
radiation ones, and y and q replaced by new ones defined as 

y = —^ ^ = ^^77 — • (44) 



eq^eq 



Prom (40), one can proceed to calculate as 

^ 3$ i S 

C = — = — H (45) 

where x — qy is <^ 1 when a mode is outside the horizon (note that it is no longer H, 
but rather Ha^ — Heqa^^, that remains constant in the radiation era). Now it has been 
demonstrated (Mahk et al (2003)) that for reheating at a constant rate of energy transfer 
from inflaton to radiation there exists a mode of constant This mode matches the C 
mode of (43) in both ( and ( at the start of reheating (the D mode is already redshifted to 
an unobservable amplitude by then) and with the A mode of (45) at the end of reheating 
provided: 

A = -C. (46) 
Consequently, for the ensuing radiation era solution one may set 

C=-p (47) 



and hence 



One then finds 



^ = -2^(^-cosa;). (48) 



2C 

5 = -2{x^' + $ + a;^$) = — (-2sinx + 2xcosx + 2x^sinx - x^cosx), (49) 
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where now = d^/dx. Thus we see that the amphtude of the fluctuations after re-entry 
(when X = l/\/3) is approximately the same as it was at horizon exit. Far ahead of re-entry 
{x <^ 1), however, the density contrast and the velocity divergence of (39) are given by 

^ = 5r = -, and Or = ^ — = - ^ C = j^C. 50 

For completeness, it can also be shown that the (unobservable) B solution of (45) and 
(29) yields the variables 

, , 2B 2B , ^ mB 3HB , , 



7. The synchronous gauge 



There are three evolution equations in this gauge that can be derived using Newtonian 
physics and the Copernican Principle (see Appendix A): 

= -19^ - Ih- 9^ + H9^ = —6^; h + 2Hh = -16nGpr6r. (52) 

The script 5* is used to distinguish the radiation density contrast 6^ and divergence of the 
peculiar velocity 9^ = ih ■ from their corresponding values in the longitudinal gauge, Sr 
and 9r, which are unscripted. The metric perturbation r] of (3) is fixed once 6^ , 9^ , and h 
are, by the relations 

4nGa^pSs Ha? ■ . _ AnGa^ 
"P ^ 2F ' ~ ~P 
which in this instance reduces to 



r] = u + -T^k V = —r^{p + p)0, (53) 



In general, a third order ordinary differential equation results from the elimination from (52) 
of two of the three perturbation variables 5^, /i, and 9^ . 

In the radiation era preceding inflation, when the modes are subhorizon, it is usual to 
assume that 

~ -^^r, 01 > \hl (55) 

in which case one can combine (52a) with (52b) to obtain the second order equation (Peebles 
(1980)) 

+ = (56) 
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The solution of (52a) and (52b) is 

5f ^= f -^dt + ^o, (57) 

J \/Sa 

and 9^ — — 3(5f/4. To solve for h, substitute (57) into (52c) to convert the latter to 



with the solution 



cth|Ei(i¥') , ho 



h^3\C\^^ + f, (59) 



where Ei[x) is the exponential integral function. 

For subhorizon modes with k ^ Ha, hence (/? ^ 1, Ei{i(fi) approaches the limit iir, so 

that 

h = ^ ^ 6inH\Sth\ + hoH. (60) 

i 

As a self-consistency check, it is easy to verify that (60) and (57) satisfy (55) provided 

l^ol « I^^^'^It^. (61) 

(here the reader is reminded of the fact that ti marks the end of radiation domination and the 
onset of inflation). Unless the universe was born with a large \ho\, however, it is reasonable 
to assume that 

\h\ ^ HlSf^l (62) 
during the radiation era between Big Bang and inflation. 

At times t S> when inflation takes over and H is very approximately constant, the 
right side of (52c) becomes neghgibly small. The solution of (52) is, for Sr and h, 

h--,+c; ^r-fcos^-^^+gj (63) 

To maintain continuity with the thermal oscillations of the radiation era of t ^ tj, it is clear 
that / ~ \6f^\ and b must be chosen to secure a smooth connection with (62) at t = tj, viz. 

\b\^aM% (64) 

which renders the last term of (63) negligible. When t > t^, (55) continues to hold, because 6^ 
for both independent modes of perturbation oscillate with the amplitude \S^\ until horizon 
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exit and becoming either a constant or fading away as a ^ afterwards (see below), while h 
is damped more quickly as a~^. 

Let us examine more closely the behavior after horizon exit. The two possible orthogonal 
solutions from (63), viz. g — or g — — 7r/2 then become, to a high degree of accuracy for 
5^ and 6^, 

(S - i^^'i - (4) i^^'i^ (I) - (4.) - (zf ) lei. (05) 

for t ^ ti- By means of (54), (64), and (65), one can estimate or calculate the remaining 
variables, viz. 

\h\ ^ = H (^Y l^^^l « H\6l!^\- \r^\ ^ = ^^lei « l^^'l, (66) 

for both pairs of (5^, 9^) in (65), and 

for the first and second pair respectively. It is then clear that 

\r}\ < ^I'^^'^l (68) 

for both pairs. The consequences of (66) through (68) will become apparent as we connect 
the two solutions here with those of the longitudinal gauge. 



8. Proof of gauge consistency 

During the inflation era when a mode is expelled from the horizon, there could in 
principle be significant differences in the perturbation variables as compared between the 
synchronous and longitudinal gauge. To check for consistency, one appeals to the gauge 
transformation equations (Ma & Bertschinger (1995)), viz. 



and 



C = V- \5l (70) 



Comparing the 5^ part of (65) with (38) one finds that 

5r = and Or = (71) 
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is secured by the simple correspondence 

C = IS'J^l and D = -ISf^l, (72) 

consistent with (35). This is because the h and fj contributions to (69) are insignificant for 
both mode solutions {i.e. both {5r, 9r) pairs), as can be verified with the help of (66) through 
(68). Moreover, with the same correspondence between C, D, and \Sl^\, (70) yields the same 
result for the curvature perturbation ( as (43), due to the smallness of rj as given by (66). 

However, the situation is a little different in the radiation era after reheating, as rj plays 
a key role in the gauge reconciliation. Here, the two solutions to (52) involve 5s — ^ 
and a are better known (see e.g. the last two of (9.121) of Kolb & Turner (1990). More 
precisely the former is 

{6^, el h,h)=Q ^ff, -f, -2i/eq) , (73) 

where Q is a normalization factor and y and q are as defined in (44). Note that the sharp 
rise of 9^ w.r.t. q (or k) explains why (52a) is dominated by the h term on its right side 
before re-entry, but by the 9^ term after. The variables rj and rj of (53) are 

One may now apply (69) and (74) to transform 6s and 9^ to the longitudinal gauge. The 
results are in agreement with (50) provided Q is set at 

Such a choice of Q would also enable one to calculate ( using the synchronous gauge variables 
and (70). This leads to C = —C/A^, the same as (47). Moreover, by (72), C, has the amplitude 
of the density fiuctuations in the thermal radiation on the horizon, at the time of horizon 
exit of the scale in question. 

The second solution 5^ = 5f ~ a is 

(5,^ , h:h) = P 2yH, -^y, -^yij) , (76) 

for some P. Consistency with the longitudinal gauge solutions for 5^, ^r, and viz. (51) 
and the B term of (45), is secured by the transformations equations of (69) and (70), after 
setting P = -3Bq/S. 
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To the best of the author's knowledge this is the first time the longitudinal and syn- 
chronous gauge are formally reconciled in detail by explicit gauge transformation. Previous 
efforts include Ma & Bertschinger (1995) who demonstrated consistency of results upon 
horizon re-entry by numerical technique, and Weinberg (2003) who was only concerned 
with the adiabatic (constant) mode at times when it is outside the horizon. 

9. Validity of lineEir evolution 

The key message of the past 3 sections is the conservation of the longitudinal gauge 
density contrast 5.^ of relic radiation, and of C) when the scale of interest is outside the 
horizon, i.e. (50), (72), and (40). Thus, if there was a pre- inflationary radiative FRW stage, 
there will be 0(1) fluctuations on scales comparable with the peak wavelength of the black 
body spectrum, i.e. /c ~ T^. By (6) one gets 

l/A;~3e^/ro. (77) 

Thus for = 60 e-folds^ this scale would eventually be inflated to a few Mpc. Long before 
that, at the time when the scale exits the horizon, its 5^ ~ 1 thermal amplitude 'freezes out' 
at this large value, when the method we used thus far involves only the lowest order of small 
quantities. 

It might be thought sufficient to reinstate linear theory that the gravitational potential 
$ Ga^hpjl? <S 1. Beware however that the sound speed = 5P/5p is an indispensable 
coefficient of (23), i.e. without precise knowledge of (23) is useless; yet the classical value 
of becomes very questionable on sub- wavelength scales when S — 5p/ p > 1 and likewise 
for P. A similar argument can be made about the growth equations in the synchronous 
gauge, e.g. (52) acquires an extra non-linear term of the form hS^ that must be included 
when > 1. The situation is worse than even this. When the relative density contrast Si of 
each fluid component is of order 1 or larger, then despite $ ^ 1 the vector and tensor modes 
arc coupled to the scalar modes and become important (Bertschinger (1993); Mollcrach & 
Matarrese (1997)). Moreover, one should fully expect very significant quantum corrections 
on the sub-wavelength behavior of fluids , as is already indicated by the change in form 
of the power spectrum on the scale of a wavelength, viz. (17) versus (19). Although it is 
conceivable that quantum and nonlinear effects might serve to smooth out the perturbations. 



^For more e- folds of inflation the observed cosmological scales would have originated from su6- wavelength 
ones, k > T, about which the density contrasts are obtainable only by the quantum field calculation, resulting 
in (19). 
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it seems far more likely that they would be observationally unacceptable. At least such a 
fortuitous outcome needs to be demonstrated, but it has not). 

For N — 70, the scale would even be larger than the Hubble radius. In that case all 
cosmological scales would be in the region where A; ^ 3], and in which our calculation gives 
perturbations with amplitude larger than unity. Exactly how this should be interpreted is 

not clear, but in any event it strongly suggests perturbations at least of order unity, again 
observationally excluded. One could appeal to the microscopic (sub-wavelength) causal 
interactions within the radiation fluid to maintain the very small scale fluctuations at their 
thermal values ahead of horizon exit. This is plausible, but is not an experimentally verified 
fact because no one has ever inflated a photon. 

Suppose one does proceed with such a classical thermodynamic argument regardless 
of initial length scale, by starting with the assumption that as the thermal wavelength 
scale exits the horizon the pre-inflationary radiation has already been redshifted to oblivion, 
Pr ^ Pv There will then be no major impact of 5r upon the overall curvature parameter C, 
if one adopts the standard 1 -\- — e ^ 10~^ equation of state for the slow-rolling inflaton 
rather than 1 -t-u;^ = 0. Then, the contribution of 5r to the curvature perturbation C, would, 
from (40), be redshifted away as in the course of inflation, and would soon become 
completely negligible {i.e. so long as C, does not become large there will not be an issue with 
the perturbation approach). To elaborate further, the relatively large magnitude oi e ^ 10^^ 
means the quantity 'iHbpIp^ dominates the 'iHbpjp contribution to in (40) because the 
redshift has rendered p^ ^ p,., thereby preventing ?)H5p/py from being blown up by too 
small a denominator. However, since 5r for a photon fluid that does not interact with others 
is conserved while the relevant scale is outside the horizon [i.e. for an isolated fluid (50) 
holds to a good approximation even when e is finite, see Mahk et al (2003)), beware that 
by the time of the coherent oscillation of the reheating process at the end of infiation there 
exists two moments per cycle when e = vanishes completely. 

To be very clear about the behavior of e at late stages, observe that during reheating 
the infiaton scalar field rolls down a steep slope of the V{(t)) versus graph towards the 
bottom of the potential well T^(0), where it is trapped and Wy undergoes rapid oscillations 
at the frequency ~ l^"(0min)- Specifically the two pressure and density of the scalar field are 
given to lowest order by 

p, = 02-V^(0), p, = 02 + i^(0); (78) 

so that Wy — 1 precisely at the potential minimum where F(0inin) = while — —1 the 
two topmost points of each cycle where = 0. In the standard treatment of the problem, this 
oscillation is regarded as a 'massive condensate' of average (w^,) = 0, i.e. cold non-relativistic 
particles capable of decaying into radiation at some constant rate T that becomes substantial 
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only when the inflaton is trapped at the potential minimum, a phase representing the final 
leg of infiation known as reheating (Bardeen & Steinhardt (1983)). 

Now the procedure of averaging whilst insisting upon the conservation of ( through- 
out the process is valid only if the oscillation of does not cause C to become large at certain 
times. With the relic radiation, however, this is no longer true. As already explained, once 
= —1 (or e = 0) the 3HSpr/pr = —35^/4 contribution to ( in (40) will be 'reclaimed' 
by the radiation, and since by (50) Sr is conserved all the while, perturbation does not hold 
for the comoving scale of one thermal wavelength that exited the horizon early, for which 
Sr ^ 1- In other words, linear growth theory cannot be deemed viable if on this important 
length scale there are times during inflation at which ( is affected by a non-perturbatively 
large term that also dominates all other terms. The fact that ( becomes large when Wy — —1 
can also be understood by writing ( as 

C = -Cv + -Cr, (79) 
P P 

and noting that and Cr are both conserved in the absence of mutual interaction between 
the inflaton and the relic radiation (and if the inflaton vacuum is mildly coupled to some 
material fluid the conservation law will still approximately be true, Malik et al (2003)) 
and while the scale is superhorizon. Thus, provided 1 + = e is flnite the smallness of 
the coefficient of ensures a negligible contribution of Cr to C throughout inflation and 
reheating; i.e. it is C ~ 1 and not Cr ~ 1 (the latter because of (40) and Sr ~ 1) that leads 
to perturbation breakdown, and this will most certainly happen ii 1 + Wy — 0. 



10. Conclusion 

A comprehensive calculation of the evolution of pre-inflationary perturbations from sub- 
horizon to superhorizon scales in the longitudinal and synchronous gauge is presented, with 
gauge reconciliation at a level of detail not easily accessible in the published literature (addi- 
tionally an Appendix is also provided to show how the growth equations in the synchronous 
gauge can be derived without appealing to the Einstein Field Equations) . It is concluded that 
the standard 'minimal' version of inflation and reheating cannot involve an earlier radiation 
dominated FRW phase. 

One obvious 'way out' is to modify standard inflation slightly. If the inflaton has been 
dissipating into another population of photons or relativistic particles even during the slow- 
rolP phase {e.g. warm inflation, Berera (1996a,b)), the 3HSp/p term of (40) would, whenever 



^Dissipation in the reheating phase itself can also serve the purpose provided a sufficient amount of fresh 
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Wy — —1, have as its denominator p — —AHpr where pr is now dominated by the much larger 
density of the decayed particles rather than the relic radiation. Thus, so long as there are 
no 0(1) fluctuations in the density of these particles among any of the scales that exited the 
horizon (which is a reasonable expectation because the new population docs not supercool, 
i.e. the large fluctuations are associated with smaller scales that do not have to exit), linear 
theory may still be robust. 

But upon further reflection this scheme docs not work, because during coherent oscil- 
lation varies between -1 and 1, so that as the scalar fleld rolls back down the potential 
well and 1 + increases away from zero the —4Hpr term will soon be offset, i.e. there will 
always be some epoch of the oscillation at which p vanishes and the contribution of the relic 
radiation to C becomes large. 

There are other ways in which the theory could adapt itself to avoid the problem. 
The universe could be born inflating (Vilenkin (1982); Hartle & Hawking (1983)). Later 
versions of inflation, such as chaotic inflation (Linde (1983, 1986)) do not necessarily begin 
with a pre-inflationary FRW phase, and would therefore be unaffected. Yet the point is 
that, as explained in the introduction, if the smoothest possible component of massless 
particles already caused such severe problems, it is hard to imagine how the presence of any 
equipartition fleld of other matter-like particles, be they curvaton or relics, would be more 
accommodating to the theory. Thus the claim that inflation also simultaneously solved such 
other major problems of cosmology as flatness and the absence of monopoles is overstated. 
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A. Newtonian derivation of the growth equations synchronous gauge 

We restrict our treatment to a radiation dominated universe. The continuity equation 
of mass and momentum conservation is 

| + V-[(p + p)v] = 0, (Al) 

where V = Vr = d/dH with 

R = ar (A2) 
as the physical distance and d/dt is done with R held fixed, i.e. 

dp 



dp 
dt 



R dt 



. dp dp 
J. Or dt 



-Hr Vr p, (A3) 



where use was made of Vr = Vr/a, and r|R = — dR/a^ = —Hr. Hence (Al) becomes 



dp 
dt 



-HvVrP+ ^ Vr • {Har + u) + ^ " • Vr(p + p) = 0, (A4) 
a a 



where we wrote 

V = R = Har + u (A5) 
with u = ar being the peculiar velocity. 

It is now necessary to discuss the diverging r ■ Vr p term, which is to do with the 
fact that in terms of the cosmic time t and the physical coordinate R = ar the metric is 
ds'^ = dt^ — {dR — HRdt)^, and has cross terms, i.e. Qqq = 1 — H^R^ and g^j = HRK That 
means there should be extra terms in the energy and momentum conservation equations. 
Fortunately, near any particular location, matter should behave more or less as it does in 
fiat space. In terms of local coordinates t and R the equations will be the same as fiat space- 
time for small R. There will be differences for large R, but if one is interested in equations 
connecting the local quantities here, one can legitimately set R = (or equivalently r = 0) 
and does not need to work out the equations for other values of R separately, because the 
cosmological principle tells us that in terms of their locally defined variables, the equations 
would be the same as they are at the origin. 

Undertaking this step, and dropping henceforth the suffix r from dp/dt and V, we 
obtain 

^ + 2>H{p + p) + ^Vr • u + - • Vr p = 0. (A6) 
dt a a 

In a homogeneous background universe the V terms vanish, and so to lowest order one 

recovers the familiar equation dp/dt — —3H{p + p). The next order is, for a radiation 
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dominated universe with p + p — Ap/3, 



Sp = -ipSH - AHp5 - ^-\ph- AHp5 - UO, 



3a 



where 



ku 



is the velocity divergence (expansion). (A7) can be compressed to read 



(A7) 
(A8) 

(A9) 



which is the same as the first of (52). In general, for a medium with any equation of state, 
(A9) becomes 



5 + 3(c^ - w)H5 = -(1 + w) I ^ + ^ J , 



(AlO) 



where = dP/dp is the speed of sound, and w — p/p. 

To get the remaining equations of (52) we turn to the Euler equation of momentum 

1 



flow: 



dv 
'dt 



+ (v • Vr)v = 

R P + P 



Vr p - Vr$, 



(All) 



where $ is the total gravitational potential here. Bearing in mind again ( A2) and ( A5) , we 
repeat similar exercises as before: 



dv 
'dt 



R 



dv 

'dt 



H{r ■ Vr)v; (v • Vr)v = H{r ■ Vr)v + 



U - Vr 



V. 



Dropping once more the suffix r from dp/dt and V, 

u- V 



(9v 



Next, enlist (A5) another time to write 



-^—Vp - -V$. 
p + p a 



and moreover 



dv 

— = ar + u + Hu] 
dt 



uV uV,^^ , uV 
V = [Har) H u. 



a a a 

The last term is second order of small quantities. The Har term gives 



(A12) 
(A13) 

(A14) 
(A15) 



u- V u- V 

V = 



H 



d d d 
(Har) — — ( ax— — h ay— — h az— 1 (axi + ay\ + az]s.) — Har — Hu. 
a a \ dx dy dz ' 



(A16) 
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(A13) now reads 

ar + u + 2Hu = —Vp - V$. (A17) 

P + P 

The background (lowest order) terms give the usual deceleration parameter q, with V$ 
interpreted as a force that depends on the nature of the cosmic substratum. The first order 
term is 

1 P 

u + 2Hu = rV(M; ot9 + H9 = —6, (A18) 

a[p + p) 4a^ 

where the last equation is the same as (52b) for any radiation era with p/p — dP/dp =1/3, 
and is obtained by means of (A8) and the synchronous gauge condition (sometimes called 
'free fall' condition) — k'^S^^ / {2Ha?). For a general medium this equation is 

d+{2- Zw)Hd = 77^^^. (A19) 
(1 + w)a^ 

Finally, (52c) is just the perturbed version of the Friedmann equation, which for the 
radiation era is 

H + H'^-^p. (A20) 

Specifically by defining SH — h/6 and taking account of the rate equation p — —AHp, one 
obtains (52c), the general form of which is 

h + 2Hh = -3H^{1 + 3cl)6. (A21) 



